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We consider the one-dimensional partially asymmetric zero range process where the hopping rates 
as well as the easy direction of hopping are random variables. For this type of disorder there is a 
condensation phenomena in the thermodynamic limit: the particles typically occupy one single site 
and the fraction of particles outside the condensate is vanishing. We use extreme value statistics 
and an asymptotically exact strong disorder renormalization group method to explore the properties 
of the steady state. In a finite system of L sites the current vanishes as J ~ L~ z , where the 
dynamical exponent, z, is exactly calculated. For < z < 1 the transport is realized by N a ~ L 1_z 
active particles, which move with a constant velocity, whereas for z > 1 the transport is due to the 
anomalous diffusion of a single Brownian particle. Inactive particles are localized at a second special 
site and their number in rare realizations is macroscopic. The average density profile of inactive 
particles has a width of, £ ~ <5 -2 , in terms of the asymmetry parameter, 8. In addition to this, we 
have investigated the approach to the steady state of the system through a coarsening process and 
found that the size of the condensate grows as ul ~ i 1 /( 1+z ' for large times. For the unbiased model 
z is formally infinite and the coarsening is logarithmically slow. 



I. INTRODUCTION 

The properties of interacting many particle systems 
can be strongly affected by the presence of quenched 
disorder. This in particular true for systems of self- 
driven particlesi^, where even pointlike defects^ are 
able to change the macroscopic properties of the sys- 
tem. Often these defects cause phase separated states, 
a phenomenon, which is known from jam formation at 
bottlenecks^. 

Particular interesting features arise, if not only the am- 
plitude of the hopping rates are quenched random vari- 
ables, but the directional bias as well. Then the dynam- 
ics of the particles is governed by a complex landscape 
of energy barriers. As the escape time growth exponen- 
tially with the heights of the barriers, the largest barri- 
ers in the system determine the velocity of the particles. 
This property of the many particle system is in agreement 
with more classical problem of single particle diffusion in 
a disordered environment^, which is rather well under- 
stood and serves for numerous applications, e.g. poly- 
mer translocation through a narrow pore& or the motion 
of molecular motors on heterogeneous tracks 7 -. 

Since there is no general framework of studying 
nonequilibrium disordered systems it is of interest to in- 
vestigate specific simple models. Here we consider the 
zero range process 8 (ZRP) with quenched disorder. The 
ZRP is particularly well suited for theoretical analysis 
because the stationary weights of a given configuration 
factorize and can be exactly calculated^. But the ZRP 
is not only conceptually interesting, there are a number 
of important applications as well: it has been used in 
order to describe e.g. the formation of traffic jamais, 
the coalescence in granular systems, and the gelation in 
networks. 

The effect of quenched disorder on the properties of 



the ZRP have been investigated in one dimension in the 
totally asymmetric versionAl, i.e. when the particles hop 
in one direction with position dependent rates. In this 
case a dynamical phase transition takes place from a 
low density phase, where one observes the condensation 
of holes, to a homogeneous high density stat o 15 i 16 . At 
low densities the average speed of the particles is related 
to the lower cut-off of the effective hopping rates. The 
properties of the phase transition are determined by the 
asymptotics of the effective hopping rate distribution at 
the lower cut-off and not by its mean or variance. If the 
ZRP is partially asymmetric and the easy direction of 
hopping is random as well, one observes a strong depen- 
dence of the average speed of the particles on the system 
size rather than on the density as in case of the totally 
asymmetric modelii. Already this result illustrates the 
qualitative difference between the two different realiza- 
tion of the disorder. 

The ZRP can be exactly mapped onto a one- 
dimensional asymmetric simple exclusion processi2ii2ii4 
(ASEP) if sites are considered as particles and masses as 
hole clusters. Disorder in the ZRP is transformed into 
particle dependent hopping rates in the ASEP, which is 
generally referred to as particle- wise disorder. In the 
ASEP one can also realize disorder which depends on 
sites. The ASEP both with p article- wis o 15 i 16 and site- 
wise disorder^iLiSiiS have been extensively studied, in 
particular for systems which are totally asymmetric. 

The partially asymmetric ASEP with particle- wise dis- 
order has been recently studied by a strong disorder 
renormalization group (RG) method 17 which provides 
asymptotically exact results for large sizes. This RG 
method has originally been introduced in order to study 
random quantum spin chains 20 , but afterwards it has 
been applied for a large variety of quantumSL 2 ^ as well 
as classical 2 ^ systems, both in and out of equilibrium 24 , 
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FIG. 1: The zero range process on a one-dimensional lattice. 
The top most particle at site i hops with rate p; to site i + 1 
and with rate to site i — 1. 



for a review see^. 

In this paper we use this RG method to study the par- 
tially asymmetric ZRP with quenched disorder. With 
this type of disorder there is a condensation phenomenon: 
almost all particles occupy one single site, whereas the 
fraction of particles outside the condensate converges to 
zero in the thermodynamic limit. Here we address ques- 
tions regarding the size dependence of different quanti- 
ties, such as the stationary current, the density profile, 
the mean density in the bulk and the number of particles 
outside the condensate. We also investigate the approach 
to the steady state through a coarsening process and the 
time dependence of the mass of the condensate. 

The paper is organized as follows. In the next section 
we will introduce the disordered ZRP and discuss its re- 
lation to the ASEP with particle disorder. In Sec II VI we 
introduce the RG method which is used to calculate the 
steady state current in the system. The density profile 
as well as finite-size behavior of the bulk density is cal- 
culated in Sec0 whereas the coarsening behavior is an- 
alyzed in Sec I VII We summarize and discuss our results 
in the final section. Some details of the calculations are 
given in the Appendices. 



II. DEFINITIONS AND PRELIMINARIES 

In the ZRP particles hop from site to site on a lat- 
tice and the hop rates depend on the departure site and 
on the number of particles (mass) at that site. With 
these conditions the stationary weight of a configuration 
is given in a factorized form, which offers an opportunity 
to analyze the steady state properties exactly^. The type 
of ZRP we consider in the present work is defined on a 
one-dimensional periodic lattice with I = 1,2, ... ,L sites 
and N particles. Particles are allowed to hop to nearest 
neighbor sites. The hopping rates are quenched random 
variables, where we denote the forward hopping rate from 
site I to I + 1 by pi and the backward hopping rate from 
site I + 1 to I by qi (see Fig. ^1 . Compared to the most 
general ZRP, we do not consider mass-dependent hopping 
rates, i.e. pi and qi are valid for all occupation numbers 
m > 1. A configuration of the system is characterized 
by the distribution of masses (number of particles), n;, 



where Yli=i n i ~ N- 

The one-dimensional ZRP defined above is equivalent 
to an ASEP where I — 1, 2, . . . , L particles are placed 
on a ring with L + N sites. The Z-th particle hops to 
empty neighboring sites with a forward (backward) rate 
of qi (pi) and a configuration is defined by the number of 
empty sites, ni, behind the particle (i.e. in front of the 
particle there are holes). 

The stationary weight of a configuration of the ZRP is 
given in a product form: 



where 
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is the canonical partition sum. The factors are given by: 

fl(ni)=9?' (3) 
and the gi satisfy the equations: 

gi(pi + qi-i) = gi-m-i + gi+m ■ (4) 

These equations are identical to the stationary weights of 
a random walker with random hop rates: pi (I — > I + 1) 
and qi-i (I — > I — 1). The solution of Eq. @ is given by: 



C 

9i = — 
Pi 



1 



L-l i 

pi+j 



i=l 3=1 



(5) 



where C is a constant. From Eq. Q one obtains: 

gt-iPL-i - giqi-i = gipi - gi+iqi = const , (6) 

which is a constant of motion. Taking the constant in 
Eq. 10 to be one we get: 
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One can obtain a number of useful results for the sta- 
tionary stated. E.g. the occupation probability, pi{n{), 



that the Z-th site contains n; particles is given by: 

ZL-l,N-n t 



Pi(m) = fi(m)- 

and the particle current reads as: 



Jl,n = (p - q) = 
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These expressions are simplified if the number of par- 
ticles goes to infinity, N — > oo, whereas there is no re- 
striction on the value of L. In this limit we obtain for 
the canonical partition sum: 
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where we label the sites in a way that gL = rnax({<7/}). 
Then the stationary current reads as: 



Jl = 1/9l, N -*■ oo , 



(11) 



and the occupation probability follows a geometrical dis- 
tribution: 

p{ni) = (1 - ai)af, ai=gi/g L , N -> oo , (12) 
so that 



(ni) =ai/(l-ai), N 



(13) 



For large, but finite N and L the sum of the occu- 
pation numbers should be equal to the total number of 
particles^: 



£<*>=Ei 



i 



1=1 



I(9iJl,n) - 1 



= N 



(14) 



In this work the hop rates are independent and iden- 
tically distributed random variables taken from the dis- 
tributions, p(p)dp and Tr(q)dq, respectively which will be 
specified later. Generally we allow for the existence of 
links with pi > qi as well as links with pi < qi with finite 
probability, i.e. the easy direction of hopping is a random 
variable, too. 

We introduce a control-parameter, <5, which character- 
izes the average asymmetry between forward and back- 
ward rates: 



[mp] av - [Inq] 
var[lnp] + varflng] 



(15) 



such that for 5 > (5 < 0) the particles move on aver- 
age to the right (left). Here, and in the following [. . .] av 
denotes average over quenched disorder, whereas var(x) 
stands for the variance of x. 

One can show, that for any non-zero density, p 
X I. > 0, a Bose condensation occurs, in the sense that a 
finite fraction of the particles, {ni)/N > 0, are condensed 
at site L, and (ni)/N tends to one in typical samples 
for N — » oo. The condensation of particles can be un- 
derstood by analyzing the distribution of g^s, a s we w ill 
show in the next section (see for comparison Refs. |l6ll"7| ). 
We shall also show that the stationary current vanishes 
in the thermodynamic limit, i.e. hin^^oo Jl = 0. These 
results are in complete agreement with a recently intro- 
duced criterion^, which predicts strong phase separation 
for vanishing stationary current. 

The numerical calculations are carried out using two 
kinds of disorder distributions. First, a bimodal distri- 
bution, where PiQi — r holds for all i and 



pip) = c5(p - 1) + (1 - c)S(p - r), 



(16) 



with r > 1 and < c < 1/2. Second, we use a uniform 
distribution defined by: 



pip) = Po 1 e(p)e( Po -p), 

Tr(q) = e(q)Q(l-q), 



(17) 



where po > and Q(x) is the Heaviside function. 

In the following we use the terminology, that the model 
is asymmetric or biased, for 5 > 0, and it is unbiased for 
5 = 0. This latter model is realized if the distribution 
of the hopping rates is symmetric. For random quantum 
spin chains, which show analogous low-energy properties, 
(5 = 0, corresponds to the critical point and (a part of) the 
5 > region is the so called Griffiths phased. Since the 
same type of mechanism takes place in the two systems 
we use the terminology of Griffiths phase for the random 
ZRP with 6 > 0, too. 



III. ANALYSIS OF THE CURRENT 

In section[n]we pointed out that the stationary weights 
of the disordered ZRP are related to the statistical 
weights, gi, of a random walker in a random environ- 
ment. For the further analysis it is important to notice 
that the random variables g defined by Eq.©, are so- 
called Kesten variables^, because there exist a number 
of rigorous mathematical results concerning their distri- 
bution Pl^). In the Griffiths phase, i.e. for 6 > and in 
the thermodynamic limit, L — > oo, the limit distribution 
presents an algebraic tail: 
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where z is the positive root of the equation 



l/z 



= 1. 



In particular we obtain for small, S, that: 



(18) 
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which is divergent and independent of the actual form of 
the distributions, p(p) and ir(q). 

If we now apply the results for Kesten variables given 
above to the disordered ZRP two remarks are in order. 
First, the algebraic tail in Eq. I|18fl should be present in 
the large (but finite) L limit. Second, the gi-s in Eq. (J5J 
are not strictly independent. However, as will be shown 
later gi and gk can be considered to be uncorrelated if 
\l — k\ > £, where £ is the finite correlation length of the 
problem. Therefore gL has to be treated as the largest 
event of a distribution in Eq. i|18|l among ~ L/£ terms. 
According to extreme value statistics^ the typical value 
of gL follows from the relation: g L 1 L = 0(1), such 
that: 
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(21) 



Another result, which we obtain from extreme value 
statistics, is the typical value of the i-th largest occu- 
pation probability, g(i)/gL- 
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From these results and from Eq. JJjJ follows that the typ- 
ical value of the current is 
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(23) 



thus vanishes algebraically in the thermodynamic limit. 
Our second result concerns the typical value of parti- 
cles outside the condensate, N out , which follows from 
Eqs.lO and (JUJ) as: 
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(24) 



Therefore, for any z > the fraction of particles in the 
condensate goes to one in the limit of large system sizes 
L — > oo and constant density, p = L/N. The typical 
number of particles outside the condensate, however, be- 
haves differently for z < 1, when it is divergent, and for 
z > 1, when it tends to a finite value. We shall discuss 
this issue in more detail in SeclVl 

Next, we turn to analyze the behavior of the gi weights 
in the weakly asymmetric limit, 6 <C 1 and estimate the 
correlation length, £. For our analysis it is convenient 
to use a correspondence between the sequences of the 
hop rates and random walk paths. (A similar reasoning 
has been introduced for the random transverse-field Ising 
chain in Ref^S.) We consider the situation where the 
walker starts at i = 0, xo = and takes in its Z-th step 
an oriented distance of 6xi = lng;_i — lnp;, such that 
its position, xi, is given by: xi = V\ 6x j, see Fig. [5] 
For 5 > the motion of the walker is biased (due to the 
average slope of the energy landscape), but decorated 
by fluctuations, i.e. one observes local deviations from 
the average behavior. The typical time during which the 
walker makes a large excursion, £, against the bias follows 
from the Gaussian nature of the fluctuations and given 
by^: 

C ~ <T 2 . (25) 




FIG. 2: Mapping the configurations to random walk paths. 



can be chosen out of ~ L positions, therefore the typ- 
ical value of follows from extreme value statistics as: 
exp(-A L /£±)L = 0(1), thus A L ~ lnLtT 1 . Keeping 
in mind the definition of Sxi and putting this result into 
Eq. we obtain: IrigL ~ 6 hiL, which is compat- 
ible with the exact result in Eq. (|21|l and the small 6 
expansion of z in Eq. (|20|l . 

In the random unbiased ZRP 6 = and the corre- 
lation length is divergent. In this case the transverse 
fluctuations of the walker in the thermodynamical limit 
are unbounded. In a finite system they typically behave 
as: Al ~ L 1 / 2 . From this follows that the fluctuations in 
the particle current in a typical sample are of the form: 



| In | J L || -L 1 ' 2 , 6 = 



(26) 



We can thus conclude that in the random walk picture 
the largest local fluctuation of the energy landscape is re- 
sponsible for the small value of the particle current. The 
particles are accumulated in front of this large barrier 
and built the condensate at a given site of the system. 
At other subleading barriers there are only a few, O(l), 
particles the sum of which gives typically N out , as esti- 
mated in Eq. pi)l. Later in Sec IV B 1 1 and I V B 31 we shall 
use this random walk mapping to obtain the size of the 
cloud of particles around the condensate and the density 
profile. 



More precisely the probability of an excursion time, i, 
reads as P(i) ~ exp(— Similarly, from the Gaussian 
form of the fluctuations follows the typical transversal 
size of excursions of the walker, which scales as, £j_ ~ 
^V 2 r*j 6~ 1 . Consequently the probability of a transversal 
size of excursions, A, is given by P(A) ~ exp(— A/£j_). 

£ as defined in Eq. (|25|l is the correlation length of 
the random walk and at the same time it is the correla- 
tion length of the ZRP as discussed in the first part of 
this section. Indeed, the stationary weight, g^, is con- 
nected to the part of the landscape which starts at po- 
sition, i, and its value is dominated by the height of the 
largest excursion, x™" ax — Xi, see Fig|5| Since the land- 
scape becomes uncorrelated for distances (times), which 
are larger than £, the corresponding gi weights are un- 
correlated, too. The value of i.e. the largest weight, 
is related to the largest possible transverse fluctuation, 
A^. In a chain of length, L, this extremal position 



IV. RENORMALIZATION AND THE 
PARTICLE FLOW 

The results about the particle current presented in the 
previous section can be obtained, together with another 
results, by the application of an RG method. Here we 
first illustrate how single sites of the lattice can be deci- 
mated, afterwards the method is adopted to the random 
system and finally the properties of the fixed points, both 
for the unbiased and asymmetric (biased) models are dis- 
cussed. 



A. Exact decimation of a single site 

The factorized form of the stationary distribution in 
Eq. allows to explicitly integrate out a typical site, 
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FIG. 3: Renormalization scheme for the zero range process: 
Site i is eliminated and hopping from site i — 1 to i + 1 and 
vice versa occurs with rates p and q, respectively. 



say i ^ L. To do this we start with Eq. © and eliminate 



Pi-iPi Qi-iQi , 

9i-i ; 9i+\ = 1, I 2 ") 



Pi + qi-i Pi + qi-i 

so that the effective (two site) hop rates are identified as 



Pi-iPi . qi-iqi 

p = —— — , q = 



(28) 



Pi + qi-i Pi+qi-i 
At the same time the particle current is transformed as: 

J = J2 ZL - 1 ' N - ni - 1 p(n l ), (29) 

rii—0 ' L 

what can be obtained along the lines of Eq. ©. From 
a practical point of view it is of importance that the 
particle current remains invariant, which happens both, 
i) for N — > oo (even for finite L) and ii) if (rii) — ► 0. For 
a random system case ii) is realized for a typical site, if 
L —> oo, even if N is finite. Indeed, in this case qi — 
0(1) whereas = N — ► oo. Since a site outside the 
condensate generally contains finite number of particles 
these are all "typical" so that one expects to be able to 
repeat the transformation up to site, L. 



B. Renormalization of the random ZRP 

In the following we shall apply the RG transformation 
in case ii), i.e. for a random system in the infinite (or 
very large) lattice limit. The corresponding single par- 
ticle problem of a random walker in a random potential 
is thoroughly studied in the literature and many exact 
results have been obtained by using the strong disorder 
RG method 2-23 . The principles of the RG procedure ap- 
ply as well to the the ZRP with many particles. There 
are, however, questions which are specific in the large N 
limit, such as the properties of the condensate, the profile 
of the particles and the coarsening process towards the 
stationary state, which will be studied in the following 
sections. 

In the following we adopt the RG rules in Eq. 1|28|) 
for the random system. In the traditional application 
of the strong disorder RG method for a Brownian par- 
ticle one renormalizes the energy landscape^ in Fig. [21 
Here we apply a somewhat different approach and also 
refer to a mapping to random quantum spin chains. The 



first step is to select the site (bond) to be decimated 
out. We choose the fastest rate, Q — max({p,}, {qi}), 
which means that processes which are faster than the 
time-scale, r = I/O, are integrated out. If the distri- 
bution of disorder is sufficiently broad, then the rate O 
is much larger than the neighboring rates and the ex- 
pressions for the renormalized rates in Eq. (|28|l can be 
simplified. For example if fl = pi we have: 



P~Pi-i, 



(30) 



and similarly for f2 = qi—\, by replacing qi and p^. It 
is evident from Eq. (|30|l that the generated new rate is 
smaller than the eliminated rate f2. Repeating the dec- 
imation procedure the energy scale is gradually lowered 
and we monitor the distribution of the hop rates, P{p, f2) 
and R(q, fi). respectively. In particular we arc interested 
in the scaling properties of the transformation at the fixed 
point, which is located at fl* = 0. This latter statement 
is in accordance with the fact, that the stationary current 
vanishes in the system. 

Using the approximate decimation rules in Eq. i|3U|) 
the RG equations can be formulated in the continuum 
limit as a set of integral-differential equations, that can 
be exactly solved at the fixed point, both at S = 02i and 
for S ^ O 3 ^. These calculations are identical to that of 
the random transverse-field Ising chain (RTFIC) and we 
borrow the results obtained for the latter model. Details 
about the mapping between the two problems as well as 
the fixed-point solutions are presented in Appendix A. 



C. The unbiased ZRP 

The unbiased ZRP with 5 — corresponds to the crit- 
ical point of the RTFIC. At this point distribution of for- 
ward and backward rates in Eq. (|A5|) are identical having 
the same exponents: ro = po ~ l/ln(0/Oo)- Here, Qq is 
a reference energy scale. The appropriate scaling variable 
is given by r/ = — (In f2 — hip) / In fi = — (lnO — lng)/lnf2, 
having a distribution p(rj) = cxp(— ?/)dry, r\ > 0. The 
length scale, Lq, which is the size of the effective cluster, 
and the energy scale are related as: 



n 



6 = 



(31) 



which shows unusual, activated scaling. In the ZRP in 
the course of renormalization huge particle clusters are 
created and the distance, X, they travel is the accumu- 
lated distance covered by the original particles that form 
the cluster, X = J2k=i x k- Finally after eliminating all 
but the last site, we obtain the accumulated distance 
traveled by all the particles during time, t. In the unbi- 
ased case [(X)] av = holds and the average mean-square 
of the accumulated displacement is given by: 



(32) 



6 



in agreement with the diffusion of a Sinai walker—. Fur- 
thermore an appropriate scaling combination between 
current (measured in a specific sample) and size is ob- 
tained as: | In | Ji||L -1 / 2 , which is compatible with the 
relation in Eq. I|26|) . 



D. The asymmetric ZRP 

The ZRP with a global bias, S > (or equivalcntly 
5 < 0), corresponds to the so called disordered Grif- 
fiths phased of the RTFIC. In this case a time scale 
t ~ fl^ exists, which separates two characteristic areas 
of the renormalization process. In the initial part of the 
renormalization, for Q > fi^, both backward and forward 
rates are decimated out, until effective clusters of typi- 
cal size, £, are created, where the correlation length £ 
is defined in Eq. 125( 1. through Eq. I31|) . In the random 
walk picture in Fig. [SJthis corresponds to eliminate the 
fluctuations of the landscape and obtain a monotonically 
decreasing curve. Now continuing the decimations for 
fl < fi^, almost exclusively forward rates are eliminated 
and the ratio of the typical backward and forward rates 
tends to zero. Consequently the system is renormalized 
to a totally asymmetric ZRP in which the distribution of 
the forward rates can be calculated (see Eq. l|A5(l ). and 
follows a pure power-law (see Eq. HA5(l h 
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where z is the dynamical exponent as defined in Eq. (|19|l . 
For comparison with Eq. i|18|) one should note that for 
the totally asymmetric model the weights in Eq. are 
~ 1/pi, consequently Eqs. I|33(l and 118f) are of identical 
form. Furthermore the effective particles become indeed 
uncorrelated in a length-scale of £. 

If effective clusters of size L are created the correspond- 
ing energy scale is lowered as: 
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(34) 



Consequently the smallest effective forward hop rate is 
given by: pl ~ L~ z . This implies the same relation for 
the stationary current as given in Eq. 12311 . 

The results presented in the previous subsections 
are expected to be asymptotically exact. Indeed dur- 
ing renormalization the distributions of the hop rates 
broaden without limits both at the critical point and in 
the Griffiths phase consequently the decimation rule in 
(|30|l becomes exact in the fixed point. 

The scaling forms of the current in Eas. (|23|) and 12(j|) 
are tested by numerical calculations in Ref— for the par- 
tially asymmetric ASEP with bimodal, i.e. discrete dis- 
order in Eq. 1(16(1. Here we performed calculations on 
the ZRP using the uniform, i.e. continuous disorder in 
Eq. (|17|l . The agreement is indeed satisfactory, both for 
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FIG. 4: Distribution of the steady state current in the ran- 
dom unbiased ZRP calculated from Eqs. 10 and (flip , a) Note 
that the distribution broadens with the size, which is a clear 
indication of an infinite disorder fixed point, b) An appropri- 
ate scaling collapse is obtained using the scaling combination 
in Eq. 
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FIG. 5: Distribution of the steady state current in the ran- 
dom asymmetric ZRP, i.e. in the Griffiths phase calculated 
from Eqs.JSJ and 1111 . a) In a log-log plot the distribution is 
shifted to the left with increasing sizes, the asymptotic slope 
of the curves is 1/z and can be used to measure the dynami- 
cal exponent, b) An appropriate scaling collapse is obtained 
using the scaling combination in Eq. I|23|l and the exact result 
for z in Eq. 11911 . 



the unbiased, (Fig. 0} and for the biased (asymmetric) 
(Fig. cases. 



DISTRIBUTION OF PARTICLES 



In order to characterize the microscopic states of the 
disordered ZRP it is instructive to distinguish between 
the (til) particles in the condensate and particles outside 
the condensate, which are either considered to be active, 
i.e. they behave as a single random walker and carry 
the current of the system, or inactive particles that are 
localized outside the condensate. 
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A. Active particles 

In a system with large finite number of sites, L, the 
number of active particles is expected to scale as, N a ~ 
L a , with an exponent < a < 1. The value of a can 
be estimated by using the relation: Jj, = vN a /L, where 
v is the average velocity of a Brownian particle in that 
landscape. In a finite system v scales as: v ~ L~ u , with 
lu = max(z — 1,0). Indeed a single Brownian walker 
moves with a constant speed for z < 1 and has anoma- 
lous diffusion for z > 1. Thus the particle current in 
a finite system behaves as: Jl ~ i a_1_w . This result 
should be compared to our analysis in Eq. I|23[) , where we 
distinguish between the z > 1 and z < 1, respectively. 



_Z. Single particle transport z > 1 

For z > 1 we have w = z — 1 and the current is so small 
that it is produced by a finite number of active particles, 
i.e. N a — O(l). In this case the accumulated distance 
traveled by all particles is simply given by: 

X~t 1/z , z>l. (35) 



2. Many particle transport z < 1 

For z < 1, however, u> = and a = 1 — z, thus the 
current in the ZRP is produced by N a ~ i 1_z active 
particles. Now the active particles have a constant veloc- 
ity, thus during time, t they travel a distance of ~ t. 

Note that the number of active particles, N a , is of the 
same order of magnitude as the typical value of particles 
outside the condensate, N out , see Eq. l|2"3jl. This is due 
to the fact that in a typical sample there are only finite 
number of inactive particles, this issue is discussed in 
details in the following subsection. 



B. Inactive particles 

Inactive particles in the random ZRP are also of two 
kinds. The first kind of inactive particles is found in a 
"cloud" , which is localized next to the condensate. The 
formation of this cloud is due to an attractive property of 
the energy landscape: particles that left the condensate 
and did not travel farther than a distance £ will typically 
turn back. The second sort of inactive particles are found 
to be localized at the subleading extrema of the energy 
landscape, i.e. they are at a site I, where gj in Eq. (J5J 
takes the next-to-leading value. As shown in Eq. H22|) 
the typical value of (n?) is of the order of one. Its aver- 
age value, however, is divergent and will be determined 
below. 
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FIG. 6: Average of the occupation probability, ctj, which cor- 
responds to the typical density of the cloud. The data were 
obtained by solving Eg. (1141 numerically for L — N = 2048 
and the disorder average was performed over 500000 samples. 
We used the uniform distribution in Eg. 1171 for different val- 
ues of the asymmetry parameter, 8, and thus for the dynami- 
cal exponent, z > 1. The data are rescaled according to 1B4I 
and using the length scale in 1361 with 1191 . 



1. Typical behavior of the cloud 

Here we consider a typical sample in which the cloud 
is due to the attraction of the condensate and well sepa- 
rated from the subleading extrema of the landscape. In 
this situation the density of inactive particles is expected 
to decay exponentially, n; ~ exp(—l/l w ), where I mea- 
sures the distance from the condensate and l w is the 
typical width of the cloud. We estimate l w by the fol- 
lowing consideration. The correlated cluster at the con- 
densate has the largest size among the clusters and its 
value, £l, follows from extreme value statistics. Using 
the distribution function of cluster sizes below Eq. i|25|) 
we obtain for its typical value: £l ~ £lnL. At this 
distance, / = £l, the typical value of the weight is 
gi = 0(1), and with Eqs. (fT3"| . JHJ), and (JHU we obtain 
ni ~ L~ z ~ exp(— £z,/Z M ). Consequently 

lw ~ - = £typ ~ ^ , (36) 

where we have made use of the scaling relationS^: 
?/Ctyp = z j where £t yp is the typical correlation length. 
The small 6 behavior in the last equation follows directly 
from Eas. (|25() and (|20|l . We have obtained thus the result 
that the typical width of the cloud of inactive particles is 
measured by the typical correlation length of the biased 
Sinai walker. 

The typical density profile of inactive particles has 
the same scaling behavior as the typical and the aver- 
age value of the occupation number, a*, as defined in 
Eq. (|12fl . For this latter quantity we have checked nu- 
merically the scaling form in Eq. for different values 
of 5 and satisfactory agreement is found, see Fig. [S] 
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2. Average number of inactive particles 

The typical value of inactive particles is finite, however 
in rare realizations it is possible to find even a macro- 
scopic number of inactive particles. This is due to the 
fact that the leading value of the weights in Eq. J5J, g^, 
and the subleading value, gj = g can be arbitrarily close 
to each other. The average value of [(nj)] av « N ia is dom- 
inated by such rare realizations, in which both {ni) and 
{nj) are macroscopic, i.e. they are of O(N). However, as 
we checked numerically, the contribution to the average 
value of such samples in which more than two sites have 
macroscopic occupation is negligible. Therefore in the 
subset of rare realizations we should consider only two 
active sites, L and I, and the unnormalized weight that 
the subleading site contains nj = n particles is in leading 
order proportional to g n g^~ n = a n g^ . The distribution 
of n is thus approximately P a (n) = (1 — a) /{I — a N+1 )a n 
for a < 1 and P a (n) = l/(N + 1), for a = l 33 . Thus wc 
obtain for the expectation value: 



(n)(a) = 



1 — a 



-(N + l) 



1 - a N + 1 ' 



a < 1 



(37) 



and (n) — N/2 for a = 1. Using the distribution 
function, p(a) we can average over the rare realizations: 

[{rij)] av = J Q (n)(a)p(a)da, which is dominated by the 
contribution as a — ► 1. Keeping in mind that the maxi- 
mal value of (n) is N/2, we can write 



Nj, 



[(nj)} 



l-2/N 



1 



-p(a)da ~ p(l)\nN. 



(38) 

provided that the distribution p(a) has a finite limiting 
value at a — > 1. As we have checked numerically this 
is indeed the case, both for the unbiased and the biased 
(asymmetric) models. Consequently the average number 
of inactive particles is logarithmically divergent, although 
its typical value is finite. Furthermore the ratio between 
the average numbers of the active and inactive particles, 
N a /Nia, tends to zero for z > 1 and tends to infinity, for 
z<l. 



3. Average density profile of inactive particles 

The average density profile of inactive particles, 
[(^/)]av, is proportional to P(l), which is the probability 
density, that the subleading, almost degenerate <? ~ #l 
is located at I = I. In the asymmetric model I can be 
either at the renormalized cluster of the condensate or at 
any other cluster. In the second case P(l) is a constant, 
which happens for I > £. This behavior is illustrated in 
the left part of Fig. □ 

In the unbiased model, as can be seen from Eq. (|25|l the 
width of the cluster diverges. Therefore I and L are al- 
ways in the same cluster, and as a consequence the prob- 
ability distribution is scale-free and a function of l/L. 




FIG. 7: Average density profiles of inactive particles for dif- 
ferent system sizes, obtained by solving Eg. 11411 numerically. 
The density was p = 1 and average was performed over 500000 
samples. Left: asymmetric model with the uniform distribu- 
tion, po = 3. The profile is constant outside the renormalized 
cluster of the condensate, the size of which is finite. Right: 
unbiased model, po = 1, rescaled according to H39H . 



The form of the probability distribution for / <C L can 
be obtained from the random walk picture in Fig. [21 The 
rare event for this process is represented by a landscape 
which is drawn by an unbiased (S = 0) random walker, 
which starts at one minimum and after I steps arrives 
to a degenerate second minimum. This means that the 
random walker is surviving (since it does not cross its 
starting position) and returns to its origin. The frac- 
tion of such random walks is given by^: P(l) ~ l~ 3 ^ 2 , 
I < L, consequently P(l) = L^ z l 2 P{l/L), where P(x) is 
a smooth scaling function, which behaves for small x as 
P(x) ~ x -3 / 2 . With this prerequisite we obtain for the 
average density profile of inactive particles for a given 
density p = N/L > 0: 



(ni)U(L) =ML)L-^ 2 P{l/L). 



(39) 



Results of numerical calculations are presented in the 
right part of Fig. \7\ which are in excellent agreement 
with the theoretical predictions. 



VI. COARSENING 

In the stationary state of the process typically almost 
all particles occupy single site. If initially the particles 
are uniformly distributed on the lattice, the system un- 
dergoes a coarsening process, meaning that the number 
of particles in the condensate and the typical size of 
empty regions is growing as time elapses. To be spe- 
cific we consider finite L and N and define the length 
scale asi£ l{t) = o~ 2 {t) = [jr J]f =1 (n|) (t)] av , which char- 
acterizes the typical number of particles at sites with a 
non-microscopic rij, or equivalently the typical distance 
between these sites- In the RG procedure the latter cor- 
responds to the typical cluster sizes (defined as the sum 
of the original links) at the energy scale f2. The growth 
rate of the typical distance, dl/dt, is proportional to the 
typical current at that energy scale. This leads to the 
differential equation: 



dZ 
dl 



(40) 
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FIG. 8: Time dependence of the coarsening length scale, /, 
in the Griffiths phase in a log-log plot. We used the bimodal 
distribution il Kill with c=0.2, r=0.5 where z = 0.5 (left) and 
c=0.3, r=0.5 where z ^ 0.818 (right), having a density p = 
3, and disorder average was performed over a few hundred 
samples. The slope of the straight line is the theoretical result 
in Eq.lgTJ: 1/(1 + z). 



At the Griffiths phase we have in (J23J, J/ ~ I z , thus the 
solution of Eq. (|4(J|1 is given by: 



(41) 



Here £ is the dynamical coarsening exponent. Note 
that C is exactly known through Eq. i|19|) and it is a con- 
tinuous function of the parameters appearing in the dis- 
tributions of hop rates. Close to the critical point £ 
becomes universal, depending only on the control pa- 
rameter: £ ps 1/(25), see Eq.(J20I). We note that the 
same result can be obtained by referring to the result 
of the RG procedure. As shown in subsection IIVDI for 
51 < the ZRP is transformed to a totally asymmetric 
ASEP with a rate distribution in Eq. (|33[) . The coarsen- 
ing of t he totally asymmetric ASEP has been analyzed 
in Refs. |llll5l34| with the result given in Eq. I|41l) . 

At the critical point £ diverges and the coarsening 
is ultra-slow and strictly at 5 = the length scale 
grows anomalously (logarithmically) with t. This type 
of growth is the so called anomalous coarsening^, which 
can be observed e.g. in spin glasses. 

The asymptotic time dependence of I can be ob- 
tained by inserting the scaling form of the typical cur- 
rent through an empty region of size I into the differential 
equation in Eq. 1)40(1. The scaling form is according i|26|) 



given by J; ~ e~ 



Eq. lj2tj|) has the asymptotic (large 



t) solution l x l 2 e cl ~ i, thus the growth of the length 
scale is logarithmically slow 



In 



lnt 



(42) 



Results of numerical simulations are presented in Fig. |S| 
for the Griffiths phase (biased model) and in Fig. [5]for 
the unbiased model. The agreement with the theoretical 
results is satisfactory. 

At the boundary of the Griffiths phase, where z = 
0, the coarsening is controlled by the behavior of the 
forward rate distribution at p — > 0. To see this, let us 
consider an arbitrary forward rate distribution for which 
p(p) ~ p- 1+1 / z o (z > 0) as p — > 0, and an arbitrary 



FIG. 9: Time dependence of the coarsening length scale, 
I, for the unbiased model using the theoretical combination 
in Eq. I)42|l . We used the bimodal distribution t\lty ) with 
c=0.5, r=0.1 (left) and c=0.5, r=0.2 (right), having a density 
p = 3, and disorder average was performed over a few hundred 
samples. 



backward rate distribution, such that S > 0. Analyzing 
HI9JI it is easy to show that ( > zq + 1 in the whole 
Griffiths phase 6 > 0, and at the boundary 



lim C(S) = z + 1, 



(43) 



which is the coarsening exponent of a totally asymmet- 
ric process obtained from the original one by setting the 
backward rates to zero. Thus, as expected, the presence 
of nonzero backward rates slows down the coarsening in 
so far as the dynamical exponent is larger in the whole 
Griffiths phase than that of the totally asymmetric pro- 
cess with zero backward rates, and the latter value is 
recovered at the boundary of the Griffiths phase. 



VII. CONCLUSION 

In this paper the prototype of a diffusive many par- 
ticle system, the zero range process, is studied in one 
dimension, in the presence of quenched disorder. This 
model system is of large importance since its station- 
ary state is in a product form therefore can be con- 
structed analytically. At the same time there exists a 
wide range of different realisations of the ZRP, and there- 
fore it can be applied to various problems of stochastic 
transport. We have analyzed the properties of the steady 
state of the disordered system and obtained many exact 
results. These results are expected to be generic for ran- 
dom, driven non-equilibrium systems and therefore can 
be useful to analyze more complicated, non-integrable 
processes, too. It also important to notice that our re- 
sults are directly related to another archetypical model of 
non-equilibrium transport, i.e. the ASEP with particle- 
wise disorder. 



A. Types of the transport 

The transport in the system is related to the form of 
the hop rate distribution and characterized by the value 
of the dynamical exponent z defined in Eq. I|19|l . We can 
distinguish four main types of transport in the random 
ZRP. 
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1. z = 

In this case the maximal value of the backward hop 
rates, q m ax, is smaller than the minimal value of the for- 
ward hop rates, p m in > Qmax- This situation is equivalent 
to the totally asymmetric ZRP, which has been previ- 
ously analyzed in the literature. In this case there is a fi- 
nite stationary current in the system, such that the num- 
ber of active particles (which carry the current) is N a ~ L 
and the particles have a finite velocity, v = 0(1). Con- 
densation of particles at a particular site is not generic, 
but appears only for hop rate distributions that are van- 
ishing fast enough at the lower cut-off c > (seei for 
details). 



2. < z < 1 

In this case q max > q m in and the current is vanishing 
in the thermodynamic limit as J ~ L~ z . The transport 
in the system is effected by N a ~ L x ~ z active particles, 
which move with a finite velocity, v = 0(1). Note that 
in the system there are infinitely many active particles, 
but their density is zero. 



3. 1 < z < oo 

The current in the system is vanishing according to 
J ~ L~ z , and the transport is effected by N a = 0(1) 
active particles, which have a vanishing velocity. Thus 
the transport is realized by the anomalous diffusion of a 
few biased random walker. 



4- Unbiased ZRP: z — » oo 

The average current is zero and the fluctuations of the 
accumulated displacement are due to a few Sinai walkers 
and given in Eq. (|32[1 . 



B. Problem with self averaging 

In the partially asymmetric random ZRP with z > 
the number of active particles, N a , and thus the trans- 
port is self averaging, i.e. the typical and the average 
values have the same scaling behavior. The inactive par- 
ticles, however, which are outside the condensate but do 
not contribute to the transport, have different proper- 
ties. Their typical number is 0(1) whereas their average 
diverges as In N. This latter value is dominated by rare 
realizations in which there is a second site with macro- 
scopic number of particles. 



C. Relation with the totally asymmetric ZRP 

The RG framework used in this paper has revealed 
a relation between the partially asymmetric ZRP with 
arbitrary type of initial disorder and the totally asym- 
metric ZRP with a hop rate distribution, which has a 
vanishing power-law tail. Indeed, during renormaliza- 
tion if the typical size of the renormalized new clusters 
becomes larger than the correlation length the new effec- 
tive particles perform a totally asymmetric motion with a 
power-law distribution of the hop rates given in Eq. 13311 . 
The exponent of this distribution is related to the form of 
the distribution in the original (i.e. partially asymmetric) 
model, see in Eq. Q19[l. In this way the problem which has 
been throughly studied for the random totally asymmet- 
ric ZRP appears naturally as the fixed point problem of 
the partially asymmetric ZRP. The power-law exponent 
then appears in a self-organized fashion. In the termi- 
nology of random systems both problems have the same 
type of strong disorder fixed point. We note that a simi- 
lar relation is encountered between the biased Sinai walk 
and the directed trap model 2 . On the other hand for the 
unbiased ZRP the singular behavior is governed by a so 
called infinite disorder fixed point. 

D. Condensation transition for limited number of 
particles 

Throughout the analysis of the partially asymmetric 
model we considered exclusively the case of a finite den- 
sity, p = N/L > 0. In this case there is always a conden- 
sate present, where one typically finds a finite fraction of 
the particles. However, if the number of particles scales 
a,s N ^ L u , < u> < 1, one expects that the condensate 
will disappear for sufficiently small values of to. Indeed, 
the transport in the system can involve N a ~ L 1-2 par- 
ticles, for < z < 1, consequently for uj < 1 — z all the 
particles carry current and the condensate is absent in 
the system. At the borderline case, with N — AL 1 ^ 2 one 
expects that the fraction of particles in the condensate 
varies with A and that there is a condensation transi- 
tion possibly at a finite value of A = A c — 0(1). First 
numerical results support this scenario of the condensa- 
tion transition, but a more detailed analysis of the phase 
transition, which is due to the slow dynamics of the pro- 
cess quite involved, is not yet completed. We are going 
to clarify this phenomena in a separate work. 

E. Higher dimensions 

The ZRP process has the same type of factorized 
steady state for any type of lattices, thus also in higher 
dimensions. Also the stationary weights involve site de- 
pendent quantities, gi, which are the stationary weights 
of the random walk on the given lattice. The renormal- 
ization described in Sec lIVI can be implemented numeri- 
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cally in this case. Similar calculations have already been 
performed for random quantum systems^. During the 
RG procedure one obtains large effective clusters hav- 
ing a complicated topology. One might ask the question 
whether the strong and infinite disorder nature of the 
fixed points as obtained in Id remains also in higher di- 
mensions. The answer to this question is negative. We 
know that the underlying random walk process has an 
upper critical dimension, d u = 2, so that the Gaussian 
nature of the random walk remains 36 (in d = 2 with log- 
arithmic corrections) even in the presence of quenched 
disorder. The same type of irrelevance of disorder is ex- 
pected to hold for the ZRP, too. 
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APPENDIX A: MAPPING TO THE RTFIC, THE 
RG EQUATIONS AND THEIR FIXED POINT 
SOLUTION 



observes a complete agreement with the correspondences: 
Pi J% and Qi <-> hi. Now properties of the RG flow for 
the ZRP can be obtained from the equivalent expressions 
for the RTFIC. 

For example the probability distribution of the back- 
ward, R(q, SI), and forward, P(p, SI), hop rates satisfy the 
set of integral-differential equations, 



dR(q,Sl) 



an 



R(q,Sl)[P(Sl,Sl) -R(Sl,Sl)] (A3) 

p(si,si) f n d q 'R( q ',n)R(^,n)^ 
J q q q 

p(p, si) [r(ci, si) - p(si, si)] (A4) 
r(si,si) [ dp'p( P ,si)p(^,si)^, 

Jn P Q 



having a solution at the fixed point: 



p q (p,si)-- 
R (q,si) 



Po (Q) ( SI 



l-po(fi) 



si \ P 

r Q (Sl) /SlV'^ 



SI \q 



(A5) 



with < p, q < SI. The value of the exponents, po(Sl) and 
ro(Sl), depend on original distributions and therefore on 
the value of the control parameter, 5. The specific values 
are given in Sec II VI 



Renormalization of the one-dimensional ZRP with 
quenched disorder as described in Sec. IIVI is equivalent 
to the same procedure of a random transverse-field Ising 
chain (RTFIC) described by the Hamiltonian: 



H 



i+i 



L 

E 

1=1 



ha 



i > 



(Al) 



where a^' z are Pauli spin operators at site I, and the cou- 
plings, J\, and the transverse fields, hi, are quenched ran- 
dom variables. During renormalization the largest term 
in the Hamiltonian, SI — max({ Ji}, {hi}) is successively 
eliminated. For example, if the largest term is a coupling, 
= Ji, the connected two sites, (i, i + 1), flip coherently 
and form an effective two-site cluster. The renormalized 
value of the transverse field acting on the spin cluster is 
given by a second-order perturbation calculation as: 



hih 



i+l 



n 



(A2) 



Decimating a strong transverse field, SI = hi, will result 
in decimating out site, i, and generating a new coupling, 
J, between the remaining sites, i — 1 and i + l. The value 
of J follows from self-duality, in Eq. (|A2(I one should 
replace hi by Ji. Comparing the decimation rules for the 
ZRP in Eq. ^ to that for the RTFIC in Eq. one 



APPENDIX B: AVERAGE OF THE 
OCCUPATION PROBABILITY 

Here we start with the asymptotic distribution of the 
(uncorrelated) Qi weights in Eq. (|18ll with the condition 
that the largest value is fixed, = G, so that 

Pi9\G)dg = ^J^e(G - g). (Bl) 

From this we obtain for the distribution of the occupation 
probabilities, ai = gi/gL, for large L as 



p(a)da ~ 

z a 



1 



-La L 



La 1 /* 



1 



1 



La 1 / 2 



The average of a is therefore given by: 
ap{a)da 

o 



da. 
(B2) 









r 


- - e -y (l - 


-I)} 


la 


.y \ 





y 



-l+z 



dy .(B3) 



In leading order of L the second term in the bracket can 
be neglected and we arrive at: 



L z ^const 



L z 



z < 1, 
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a 



L 1 (const + In L) ~ L l hiL z = l, independent of 5. 



[a] av « L"« (^const + -_ j ~ -_ z > l.(B4) 

Thus, as far as z > 1 the occupation probability is non- 
self-averaging, and the decay exponent is becoming one 
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